Marginal pinning of vortices at high temperature 
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We analyze the competition between thermal fluctuations and pinning of vortices in bulk type II 
superconductors subject to point-like disorder and derive an expression for the temperature depen- 
dence of the pinning length L C (T) which separates different types of single vortex wandering. Given 
a disorder potential with a basic scale £ and a correlator Kq(u) ~ Ko (u/£)~^ ln a (u/£) we determine 
the dependence of L C (T) on the correlator range: correlators with /3 > 2 (short-range) and /3 < 2 
(long-range) lead to the known results L C (T) ~ L c (0) exp[GT 3 ] and L C (T) ~ L c (0)(CT) (4 +' 3)/(2 -' 3) , 
respectively. Using functional renormalization group we show that for f3 = 2 the result takes the 
interpolating form L C {T) ~ L c (0) exp[CT s ^ 2+a ^]. Pinning of vortices in bulk type II superconduc- 
tors involves a long-range correlator with fl — 2, a = I on intermediate scales £ < u < A, with £ 
and A the coherence length and London penetration depth, hence L C (T) ~ L c (0) exp[CT]; at large 
distances L C (T) crosses over to the usual short-range behavior. 

PACS numbers: 05.20.-y, 64.60. Cn, 74.60.Ge 



Many properties of type II superconductors derive 
from the interaction between vortices and pinning 
centersBEI. Impurities are particularly important in situa- 
tions where a transport current is applied as they provide 
the necessary pinning force compensating the Lorentz 
force acting on vorticesa and thus enable a dissipation- 
free current flow. In this note we address specific aspects 
of the high temperature pinning behavior of vortices sub- 
ject to weak point-like disorder. The results are of par- 
ticular relevance for the copper-oxide high T c supercon- 
ductors which can be operated at high temperatures and 
where such a pinning landscape naturally derives from 
oxygen vacancies. 

A single flux line in apdisorder potential is described 
by the partition functionQ 



Z(u,L)= / V[u'(z)} 

exp{-I^[|(^)%V>'(zM 



(1) 



with (c/2)(<9 z ii) 2 the elastic energy and V(u, z) the dis- 
order potential which might take negative values. The 
disorder is chosen to be a Gaussian random variable with 
zero mean and a correlator 



(V(u, z)V(0,0)} = K {u)6(z), 



(2) 



where (. . .) denotes the average over disorder realizations. 
The correlator Ko(u) decays on a length £; in most appli- 
cations the function Kq(u) is assumed to be rapidly de- 
caying, and at sufficiently high temperatures the physical 
behavior is determined by the integral A = J d n uKo(u) 
alone. In this paper we draw attention to the situation in 
disordered type II superconductors where the correlator 
Kq(u) describing the potential landscape of vortices is 
long-ranged, 



P u 

K Q (u) ~ K ^ln-, 



(3) 



in the intermediate asymptotic regime £ < u < A (here, £ 
and A denote the coherence length and the London pen- 
etration depth of the superconductor). The long-range 
tail (||) renders the correlator non-integrable and modifies 
the pinning characteristics as compared to a short-range 
correlated disorder landscape. This effect is particularly 
pronounced for vortices in bulk high-T c superconductors 
where the disorder becomes marginal at high tempera- 
tures and the ratio A/£ is large (below we describe the 
situation in an isotropic material; effects of anisptropy 
can be accounted for within the scaling approachB). 

The following discussion is not restricted to single vor- 
tex pinning; rather we consider the latter as a specific 
realization of the directed elastic string (or polymer-) 
problemQ which describes such diverse physical systems 
as domain walls, in magnetic filmstn3 n wetting (in the 
plane)u, vorticesB or random polymersO. The numerous 
non-trivial features that these systems have in common 
derive from an intricate interplay of elasticity, disorder 
and thermal fluctuations. While the elastic forces tend 
to stretch the string, the disorder potential favors config- 
urations deviating from a straight line in order to take ad- 
vantage of the potential valleys. Within the weak collec- 
tive pinning scenarioli!3 the elasticity dominates on scales 
smaller than the cross-over scale L c and forces the string 
to stay in the same valley, whereas on larger scales the 
string effectively divides up into segments of size L c which 
adjust independently to the disorder landscape. Increas- 
ing the disorder strength decreases L c , while thermal fluc- 
tuations tend to smooth the disorder landscape implying 
an increase of L C (T) with temperature. Also, the typical 
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barriers separating adjacent valleys are reduced by ther- 
mal fluctuations leading to a peculiar form-of the creep- 
type dynamics under a small external forcetil. Below we 
first discuss some general properties of vortex pinning in 
disordered type II superconductors and derive the asymp- 
totic form (J3|) of the correlator. Second, we calculate the 
pinning length L C (T) using—the functional rcnormaliza- 
tion group (FRG) approachE3. 

The behavior of random directed polymers strongly 
depends on the number n of transverse motional de- 
grees of freedom. E.g., the (low temperature) rough- 
ness as characterized by the wandering exponent = 
limi^oo <91n(it 2 (L)}/<9mL 2 decreases with n (here, L de- 
notes the length of the polymer segment and u(L) is the 
relative transverse displacement of its end points). Upon 
increasing the temperature a phase transition (roughen- 

f transition) is known to occur for n > 3, see Ref. 
the high temperature phase being dominated by 
thermal fluctuations on all length scales with a thermal 
roughness exponent £th = 1/2. On the other hand, for 
n = 1, 2 the large scale behavior remains dominated 
by disorder with the same exponent £„ as at low tem- 
peratures, however, going beyond the so-called depin- 
ning temperature Td p , the crossover scale L C (T) increases 
rapidly with temperature. In the physically important 
case n = 2 the crossover length L C (T) grows exponen- 
tially, L C (T) ~ L c (0)exp[C(T/T dp ) v ], reflecting the fact 
that n = 2 is the lower critical dimension of the rough- 
ening transition. While v = 3 is a well established re- 
sult describing the situation for a short-range disorder 
potentials, the exponent is modified by the long-range 
tail of the potential correlator as it appears in the vortex 
problem: we will show below that v = 1 (a similar ef- 
fect is found in the context of individual vortices pinned 
onto columnar tracksa). We emphasize that the asymp- 
totic behavior (||) does not influence the roughness of the 
polymer on large scales, that is, the value of the wander- 
ing exponent remains unchanged. While for n = 2 the 
non-integrability of the function Kq(u) with an asymp- 
totic decay slower than l/[u 2 ln(u)], is sufficient to change 
the value of the exponent v, the criterion on the asymp- 
totics K (u)~GfcJ(o{£, / to change the value of £2 has 
been arguecOlLj to be (3 < 3/^2, sr — 4. Inserting the 
numerically known value f2,sr ~ 5/8 for the short-range 
wandering exponent, we see that the (2 = C2,sr is unaf- 
fected by the weak non-integrability (/3 = 2) in the vortex 
pinning problem. 

We briefly derive the form of the long-range cor- 
relations occurring in the problera of single vortex 
pinning in type II superconductorsB. The Ginzburg- 
Landau equation for the macroscopic wave function \E' = 
\J p(R) exp(i(p) takes the form 



2m \ 2 

$0 / 



(1 - i*! 2 )* = 0, 



(4) 



asymptotic regime. We concentrate on the vortex solu- 
tion where the phase turn in tp by 2tt drives a circular 
vortex current. On scales R < A we can ignore trans- 
verse screening and using (V(p) 2 — 1/R 2 , the real part of 
(|) simplifies to (-<f/i? 2 + 1 - |*| 2 )* = 0, hence 



|1-| 2 (£ < R < A) » 1 -£ 2 /i? 2 



(5) 



the circular vortex current produces an order parame- 
ter suppression decaying only slowly at small distances 
R < X (transverse screening quenches the current flow 
beyond the screening length A and the suppression of the 
order parameter is exponentially smalll- The same re- 
sult follows from a variational AnsatzOEJ, \& = R/(R 2 + 
2£ 2 ) 1 / 2 exp(i(p). In high temperature superconductors A 
is much larger than £ and the slow decay of the order 
parameter extends over a wide region. 

To fix ideas, we consider disarder in the critical tem- 
perature (so called 5T c -disordercl) described through spa- 
tial variations 5a (r) in the Ginzburg-Landau parameter 
a with correlations (8a(r)5a(r')) = 7<5 3 (r — r). The pin- 
ning energy of a vortex aligned along the z-axis and with 
coordinates (u(z),z) is given by 



E pia (u,z) 



d 2 Rp(\u - R|)fo(R, z), 



(6) 



where A denotes the vector potential, $0 = hc/2e is the 
flux unit, and we have normalized to unity in the 



with the vortex form factor p(R) = 1 — ^(i?)) 2 = 
2( 2 /(R 2 + 2£ 2 ). Within the range £ < u < X the pinning 
energy correlator assumes the form 

K (u) ee <£ pin (u, z)E pin {0, 0)) (7) 
= 7 S(z) f d 2 Rp(\u--R\)p(R) ~K 6(z)^]n^, 

with Kq — 27r| , I , o| 4 7£ 2 (we have assumed l^ol = 1 in 
the derivation above). The integral over this correlator 
diverges logarithmically and is cut off only at the large 
scale A. This long-range feature will have an important 
effect on the cross-over scale L C (T) at high temperatures 
which we are now going to calculate. 

We analyze the system described by the partition func- 
tion ([l]) with the help ctf-the functional renormalization 
group (FRG) approach!! 2 ! Applying momentum shell 
renormalization to the replicated Hamiltonian leads to 
the following system of one loop equations (see Ref. |i~9| ) 
for the renormalizcd correlator Ki and temperature XJ, 

d l K l (\u\) = [3-(4 + n)<;]K l (\u\) (8) 
+ CV-(u^(|u|))+f^(|u|) 

+ 1 [*niui)jrnH)/2 - Kr(\u\)Krm , 

d l f l = (l-2()f l , (9) 

where / = 1/(ttc 2 A 3 ), = Ti/(ncA), and A" 1 denotes 
the short-scale cut-off of the theory. In (|§|) we limited our- 
selves to the two replica correlator AT/, neglecting higher 
replica terms generated at high temperatures during mo- 
mentum shell integration, see Ref. 19; an analysis includ- 



ing these terms reveals that their feedback to the flow of 
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Ki is of the same order as the non-linear terms in (g) 
and hence the qualitative results obtained below remain 
valid. 

At zero temperature the fourth derivative of the cor- 
relator K^{Q) diverges at the finite scale l c (T = 0) w 
m{9/[(n+8)Ii^g (0)]}, signalling the emergence of a non- 
analyticity in the correlator at the origin and indicating 
that perturbation theory in the disorder breaks down be- 
yond this scale. Still, the FRG-flow remains well-defined, 
and for the case of short-range disorder the correlator can 
be shown to flow rapidly towards a non-analytic disorder^ 
dominated fixed point beyond Z c (0). The crossover scalecl 
L c (0) = A"V<=(°) rj (c 2 /A:< 4) (0)) 1/3 is naturally inter- 
preted as the typical length of independently pinned vor- 
tex segments within the weak collective pinning theory; 
the same result is obtained from a simple scaling argu- 
ment equating the elastic and disorder energies. 

At finite temperatures the correlator remains analytic 
due to the thermal smearing introduced via the term 
TK^ . However, we can still identify the crossover scale 
l c (T) as the value of I where the non-linear terms in 
Eq. (@) become of the same order as the linear terms, 
indicating that the disorder cannot any longer be dealt 
with perturbatively; beyond this scale K\ is again driven 
towards the (T = 0) strong coupling fixed point. 

At small scales, the non-linear terms in the flow equa- 
tion (||) for the correlator K\ can be neglected. Fur- 
thermore, since we want to study high temperatures, 
we choose the roughness exponent to take the thermal 
value C = Cth = 1/2 which is convenient since the physics 
appears more transparently in the sequel, in particular, 
the temperature does not renormalize when £ = Cth and 
T; = T (however, note that from a mathematical point 
of view the physical results below are independent of this 
particular choice of_£). The linearized flow equation can 
be solved explicitljO with the result 



Ki(u) 



i/2)l 



</2 



d"u exp 



4ttT(1 - e~ l ) 

(u-u'e-'/ 2 ) 2 



(10) 



4T(1 - e- 



K Q {u'). 



Let us analyze the situation for short- and long-range 
correlated disorder separately. If the disorder correlator 
is short-ranged, the total weight A = J d n u' Kq(u') is 
finite and the expression ( [Tof ) can be simplified at high 
temperatures f > T dp = fe' 1 "^, 



Ki(u) = e ( 1 -"/2)i A j--n/2 



exp[ 



2 /4f], 



(11) 



for max[l,ln(C 2 /f)] < I < Z c (0). Note that for n = 1 the 
correlator grows exponentially, reflecting the relevance of 
the disorder, whereas in the marginal case n = 2 the cor- 
relator flows to a fixed point. In this case we can interpret 
the linearized flow equation as the Fokker-Planck equa- 
tion for the probability distribution Ki (u) of a particle in 
a harmonic potential V(u) = u 2 /4, where I plays the role 



of time and T is the temperature. The above fixed point 
corresponds to the stationary Gibbs-Boltzmann distribu- 
tion which is rapidly approached at high temperatures. 

In the case of long-range correlations the integration 
kernel in Eq. (|l(]) cuts the non-integrable tails of the cor- 
relator at u' ~ (e'T) 1 / 2 and we may estimate the renor- 
malized correlator at the origin as 



Ki(0) 



,(l-n/2)l 



(4T) 



-n/2 



(e'T) 1 



Koiu)^- 1 du. (12) 



In the case n = 2 this expression can be understood in 
terms of the diffusive motion of "particles" that are ini- 
tially distributed with a density proportional to Kq[u). 
Their total number is infinite due to the non-integrable 
tails of Kq(u) (sec Eq. (0)), and a growing number of 
"particles" from further and further away will eventually 
accumulate at the origin, continuously increasing Ki(0) 
(see Fig. |l|) . Thermal equilibrium is never reached in this 
situation. 
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FIG. 1. The high temperature (f = £ 2 ) diffu- 

sive flow of the correlator Ki(u) evolving under the lin- 
earized FRG for n = 2. We choose a bare correlator 
K (u) = {l+ln[l+(u/£) 2 ]/2}/[l + (u/£) 2 ] with non-integrable 
logarithmic tails as in the vortex problem. After a crossover 
around / = 1 the correlator Ki (0) increases quadratically with 
I as expected from Eq. (13). The inset schematically illus- 
trates the interplay of the diffusive inflow of "particles" (hor- 
izontal arrows) and the "particle" sources originating from 
the quadratic one-loop terms (solid dots). In the short-range 
case the first mechanism is absent and Ki grows due to the 
non-linear terms only, while in the strongly long-range case 
(J3 < 2), the non- linear terms can be neglected and Ki grows 
due to the inflow from the tails. The vortex pinning problem 
is intermediate with a crossover at a scale l CI (T) < l c {T) from 
diffusion driven growth to source driven growth. 

In the following we first discuss the situation where the 
linearized flow gives rise to an exponential growth of the 
renormalized correlator at the origin. This is the case 
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for n = 1, or for n = 2 if the initial correlator decays 
asymptotically as Ko(u) ~ Ko(£/lif with an exponent 
j3 < 2. Under these conditions, the contribution of the 
non-linear terms to the growth of Ki can be neglected, 
and the crossover scale l c (T) is found from comparing 
linear and quadratic terms in the flow equation, i.e., 



T. 



(13) 



The second derivative of K\ can be estimated from 
Eq. ph, Kl'(Q) ~ Ki(0)/f. Using expression @ in the 



above crossover condition we find 



L C (T > T dp ) ~ L, 



HI) 



T \ (4+/3)/(2-/3) 
dp' 



with the depinning temperature 



dp 



io(0)' 



(14) 



(15) 



A short-range caKJ-elator with n = 1 requires /3 = 1 
and one recovers 2 !] the well-known dependence L C (T) ~ 
I/ c (0)(T/T dp ) 5 . The result (||) can also be derived from 
simple scaling arguments, c.f., Refs. p|,|il"|. 

The marginal case n = 2 is more subtle if the correla- 
tor is short-ranged or only weakly non-integrable with an 
asymptotic behavior K (u) ~ K (£/u) 2 ln a (u/£), where 
a > — 1. In the short-range case the "Boltzmannian" 
fixed point ( [Tl"| ) is unstable only due to the quadratic 
loop corrections in the FRG. They can be considered as 
sources of "particles" in the picture introduced above (see 
Fig. |). Well below the crossover l c (T), these new parti- 
cles will quickly thermalize, and the shape of the distribu- 
tion will remain Boltzmannian, whereas its total weight 
will grow. Qualitatively, we obtain l c (T) as follows: From 
Eqs. ([!(]) and ( |ll"l) Ki is seen to vary on a typical scale 
T 1 / 2 thus the quadratic terms on the right-hand side of 
Eq. (§) can be estimated by AI[Ki(0)/T} 2 for u < f 1 / 2 
with A a numerical prefactor. The amplitude of the cor- 
relator is determined from the flow equation 



1 z 



(16) 



Integrating over the interval max[l, ln(£ 2 /T)] < I < l c (T) 
and using the crossover condition ( p"3] ) as well as the ini- 
tial shape Eq. ([ll]) one recovers the well-known result 
L C (T) ~ L c (0)exp[C(T/T dp ) 3 ]. Here and in the follow- 
ing we make use of the fact that the details of the initial 
condition Kq(u) are quickly washed out, giving way to 
a correlator of Boltzmann shape (|Tl|); the latter serves 
as our new initial condition. For the short-range case an 
exact flow equation for the total weight can be written 
down under the assumption that the Boltzmannian shape 
of the correlator is preserved. This allows to determine 
the constant C — 32/tt exactly, see Ref. [l9|. 

Let us then discuss the case of a weakly long-range cor- 
relator with a tail K (u) ~ K (£_/u) 2 \n a (u/£). Note that 



now the total weight of the correlator diverges and the 
linearized flow does not approach a fixed point. However, 
the flow of the central part of the correlator is still gov- 
erned by the linear terms in the flow equation (|J) which 
force the correlator to maintain an almost Boltzmannian 
shape of width (4T) 1 / 2 . We will therefore assume that 
the flow in the region u < T 1 / 2 is captured by the growth 
of the amplitude Ki(0), and that second derivatives may 
be replaced up to numerical factors by l/T. In the first 
stage of the FRG-flow, the growth of Ki(0) is dominated 
by the diffusive flow of particles to the origin as described 
by (|l2l). However, at an intermediate scale l a (T) < l c (T) 
the quadratic source terms become non-negligible and fi- 
nally dominate, the subsequent flow being analogous to 
the short-range case. In order to determine the scale 
l rr (T) we compare the rate of growth due to diffusion, 
with the magnitude of the quadratic terms in 




5/^ diff (0) 



(12) 



[*o(i 



l M= ( e lcrT) 1 /2 



8T 



f 



(17) 



Assuming the bare correlator to decay asymptotically as 
K (u) ~ K 2 ln Q (u/£) and evaluating the right- 
hand side of ( Jl7| ) with the help of Eq. (JTJ) we find that 



UX) ~ (S-) 



T \3/(2+a) 

Xdp y 



(18) 



Note that we have to require that a > — 1, otherwise the 
inflow of "particles " will saturate already at small scales 
and the crossover condition ( p~7| ) cannot be applied. 

The further flow will be dominated by the non-linear 
terms while the linear terms still make sure that the 
overall shape of Ki resembles a Boltzmannian of width 
(4T) 1 / 2 . Only beyond the crossover scale l c (T) will the 
non-linear terms be large enough to drive the correlator 
to a disorder dominated fixed point. As in the short- 
range case, Ki(0) satisfies the approximate equation jl^ ) 
in the region l CI (T) < I < l c (T). Integrating the latter 
beyond l C r{T) and using the condition ( p^ ) we find 



Z c (f)-Z cr (f)~ (^-) 



T \3/(2+a) 



and, recalling I = ln(AZ), 



L C (T) ~ L c (0) exp 



T \3/(2+a) 



(19) 



(20) 



with an unknown numerical factor C; its determination 
as well as the precise value of the prefactor is beyond 
the accuracy of the present analysis. The short-range 
result with an exponent proportional to T 3 is recovered 
for a = — 1 which corresponds to the limiting case of an 
integrable correlator (a < —1). Eq. ( ptj| ) shows that at 
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a given (high) temperature the crossover scale L C (T) is 
the smaller the larger is the range of the potential corre- 
lations, turning finally to a power law (|lj) for strongly 
long-range correlators with f3 < 2. Note that while both 
the long- and short-range results derive from cutting the 
flow through the condition (13), the discussion of the 
intermediate range correlator with fj = 2 and a > — 1 
involves the additional crossover length l CI (T) which de- 
rives from the condition (JTt]). The latter condition then 
determines the exponent 3/(2 + a) providing the smooth 
interpolation between the long- and short-range results. 

Single vortex pinning is described by a — 1 which leads 
to an exponent oc T. However, in the above derivation 
we have not taken into account that the logarithmic tails 
(0) only extend up to A. The result ( ^o|) only applies if 
at l CT (f) the integral ( |i"o| ) is still cut by the Boltzmann 
kernel such that the estimate ([12]) remains valid. We 
may obtain an upper bound T u on T by requiring that 
( e i C r(T„)jyi/2 _ ^ j n crossover At tempera- 

tures higher than T u a crossover to the short-range case 
will take place. Finally, the range of validity of (20) is 
found to be 



T dn <T< T„. 



i dp 1 ^ J «~7d p hr» ^=j dp j 

Apart from the restriction on the temperature the mag- 
netic field B has to be weak such that effects due to 
vortex-vortex interactions can be ignored; the corre- 
sponding condition L C (T) 2 <. &o/B derives from com- 
paring tilt and shear energieso. 

In conclusion, we have determined the temperature de- 
pendence of the collective pinning length L C (T) for a di- 
rected elastic string subject to a long-range correlated 
disorder potential. In the physically important case of 
single vortex pinning in three-dimensional bulk mate- 
rial L C (T) exhibits an exponential sensitivity to temper- 
ature: in the intermediate temperature range T dp < T < 
(In k) Tdp the logarithmic tails of the correlator produce 
a simple exponential law L C (T) ~ L c (0) exp[C(T/T<j p )], 
while for T > (In k) T dp the usual short-range result 
L C (T) ~ i c (0)exp[C(T/T dp ) 3 ] holds. 



'(«) = T dp ln(«). 



(21) 
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